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PREFACE 


This report summarizes the development of theoretical rotor blade 
structural models for designs based upon composite construction. Care has been 
exercised to include a member of nondassical effects that the author's previous 
experience indicated would be potentially important to account for, A model, 
representative of the size of a main rotor blade, is analyzed in order to assess the 
importance of various influences. The findings of this model study suggest that for 
the slenderness and dosed cell construction considered, the refinements are of 
little importance and a dassical-type theory is adequate. The potential of elastic 
tailoring is dramatically demonstrated, so the generality of arbitrary ply layup in 
the cell wall is needed to exploit this opportunity. 
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INTRODUCTION 

Composite material systems are now the primary materials for helicopter 
rotor system applications. Bearingiess rotor designs proposed for the LHX 
helicopter are an example. In addition to reduced weight and increased fatigue 
lif^, these materials provide designs with fewer parts which means increased 
service life and improved maintainability. Also, in terms of manufacturing, it is 
possible to achieve more general aerodynamic shapes including flapwise variation in 
planform, section and thickness. 

The aeroelastic environment in which rotor blades operate consists of 
inertial, aerodynamic and elastic loadings. Because of the directional nature of the 
composite materials, it is possible to construct rotor blades with different ply 
orientations and hybrid combinations of materials exhibiting coupling between 
various elastic modes of deformation. For example, if the fibers are placed 
asymmetrically in the upper and lower portions of the blade, there will be a twist 
induced by flapwise bending. This provides a potential for improving the 
performance of a lifting surface through aeroelastic tailoring of the primary load- 
bearing structure. Aeroelastic tailoring of a composite structure involves a design 
process in which the materials and dimensions are selected to yield specific 
coupling characteristics which in turn enhance the overall performance of the 
structure. The design of such advanced structures requires simple and reliable 
analytical tools which can take into consideration the directional nature of these 
materials. In this report, a comprehensive bending theory is presented to aid in the 
design of composite rotor blades. 

Most of the existing analyses are formulated for isotropic metal blades. A 
general tendency is to extend these analyses to composite blades with bending and 
shear stiffnesses appropriately modified. This, however, does not describe the 
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anisotropic character of the composites. One of the earlier efforts to account for 

this particular feature of composite wings is made by Welsshaar^. He has developed 

an engineering theory incorporating the bending-torsion coupling effect for the 
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study of divergence of swept forward wings. Later Mansfield and Sobey have 
presented an account of the stiffness characteristics of a cylindrical tube 
representing a helicopter blade composed of a number of plies of arbitrary layup. 
This theory is derived within the context of Batho-Bredt engineering analysis of thin 
walled structures and small displacements. 

The above theories appear to be adequate for predicting the overall 
response of slender blades. These theories, however, do not have the capability to 
predict the stresses or the response with adequate precision when the blade is 
not slender. A complex three-dimensional stress field develops in such blades 
which is further complicated by the anisotropy of the materials and inhomogeneity 
of blade construction. Finite element analyses are reported in Refs. (3) and (4) to 
obtain stiffnesses and stresses for thick walled blades. For an accurate 
determination of the stresses, one needs to employ a large number of elements 
which makes this type of analyses computationally unacceptable for design 
purposes. 
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Hong and Chopra have studied the effects of fiber orientation on the 
flutter of rotor blades. This analysis is based on the nonlinear kinematics of 
Hodges and Dowell^. 
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Recent theoretical research has contributed a new appreciation for other 
nonclassical effects in addition to transverse shear deformation in bending- related 
behavior. A nonclassical contribution to axial stress and transverse normal strain 
affect response in planar bending situations significantly for certain combinations 
of geometry and stiffness for homogeneous structures. Incorporating these 
influences, theories were developed for homogeneous plates , stiffened plates and 
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laminates* 0 ’**. Figure 1 depicts the response contributions incorporated in these 
theories. Stresses and displacements are improved beyond linear {with respect to 
thickness coordinate) distributions, yet the overall equations retain the character 
of an engineering theory. This is a distinguishing feature for it allows the 
integration of nonclassical influences on a simple basis. 

The nonclassical influences relevant to rotor blades are those due to 
transverse shear, bending-related warping, stretching-related warping and torsion 
related warping. Laminated composites are in general strong and stiff in the plane 
of lamination and weak and flexible in the transverse direction. Consequently, 
transverse shear deformation becomes much more pronounced. Bending-related 
section warping also affects reponse in a similar way, but it is due to the fact that 
balding strain does not strictly correspond to planar deformation. Torsion-related 
warping arises whenever a section is restrained against out oi; plane deformation. 
The key to improving the stress predictive capability of a theory is to account for 
these effects correctly. With this background, the present work is undertaken. A 
comprehensive theory is developed within the context of small displacements for a 
single cell composite blade model. The distributions for the warping displacement, 
axial stress and shear flow are improved. 

The starting point is the engineering theory of Mansfield and Sobey*. The 
kinematic procedure to develop the equations follows that of Ref. (8). An example 
is presented to illustrate the potential for aeroelastic tailoring and the significance 
of nonclassical influences for rotor blades. 

OVERVIEW 

A typical composite rotor blade cross section appears in Figure 2. For a 
preliminary design, it is common to consider the forward structural box as the 
primary stiffness producing and load bearing structure. The remaining structure is 
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considered only as a mass contribution, along with tuning weight which is added to 
favorably alter the dynamic characteristics. Therefore, the foregoing development 
is for a single cell beam which is shown in Figure 3 together with a reference 
coordinate axes system. 

The classical development of the theory of thin walled beams under 
combined loading is based on two hypotheses. The first one is the Euler-Bernoulli 
hypothesis which states that plane cross sections normal to the axis of the beam 
remain plane during bending-related deformation and that these sections remain 
undistorted in their own plane. The second one is due to Saint-Venant which states 
that the applied torque is carried by a uniform shear flow and the cross sections are 
free to warp in torsion-related deformation. According to the first hypothesis, the 
only significant stress is the axial stress, a . This stress is effective in carrying 
bending and stretching-related loads. For thin walled beams, these assumptions 
lead to the following results: 


u(x, y, z) = U - yV - zW 2w <h 

fA |A jA 


( 1 ) 


where 


v(x, y, z) = V - z<|> 


w(x, y, z) = W yq> 





(la) 

(lb) 

( 2 ) 

(2a) 


Here, u, v and w denote the components of displacement of any point on the cross 
section in the directions of the coordinate axes x, y and z, respectively. Similarly, 
U, V and W represent the respective displacement components of the origin of the 
(x, y, z) coordinate system. The twist of the cross section is denoted by <j> . In 
Saint-Venant-type torsion, the variation of this variable is linear with respect to 
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the axial coordinate, x. U, V, W and $ are functions of x only. It is convenient to 
describe the variation of any given parameter on the cross section with respect to a 
circumferential coordinate, s, whose origin can be arbitrarily selected. As shown in 
Figure 3, the perpendicular from the origin of the (x, y, z) coordinate system on to 
the tangent at any point s on the circumference is denoted by r. N is the shear 
flow. It is the shear stress resultant per unit length of circumference. M x is 
the applied torque and ft is the area enclosed by the box. 

There are additional assumptions made regarding the stress state in the 
structure. Let s and n denote circumferential and tangential coordinates so that (x, 
s, n) forms an orthogonal coordinate system. With respect to this coordinate 
system, o v and a are the dominant stresses in the structure. The axial stress, 
a, is due to bending and stretching and the shear stress, a , is due to the applied 
torque. In the classical development, the remaining stresses -- a R , cr s , a ns and cr^ 
are assumed to be zero because their influence on the response is negligible. This 
is justified for the following reasons. 

In view of the nature of the construction, the structure cannot develop and 
support a ns . Therefore, this stress is neglected. The skin of the structure is 
usually very thin. Therefore, the skin is assumed to be in a state of plane stress 

with respect to the n-coordinate direction, which means that o and o are zero. 

r n nx 

It is essential for the aerodynamic reasons to preserve the airfoil shape. The 
rigidity of the cross section is maintained by an internal structure, such as a 
honeycomb filling which does not add measurably to the overall stiffness of the 
structure. Therefore, the circumferential normal stress, o ^ cannot develop. It is 
also not possible to account for this stress within the- limitations of a one- 
dimensional theory. In the new theory, these stress assumptions are retained. 

In order to account for the transverse shear, bending-related section 

warping and torsion-related section warping, a priori knowledge of a and o 

X xs 


is necessary. Statically equilibrating stresses can be developed based on the 
assumptions given by Equations '') and (2). As a first approximation, it is assumed 
that these stresses are sufficient to estimate the nonciassical effects. A brief 
summary of the procedure for developing these stresses follows. This procedure is 
similar to that of Mansfield and Sobey 1 up to Equation (16). It serves to introduce 
notation and provides a basis fa- the new theory. A flow diagram of the subsequent 
development of the new theory is shown in Figure 4. 

DEVELOPMENT OF ENGINEERING STRESSES 

The skin of the single cell tdbe is considered to be a laminate built from 
distinctly different layers or lamina. For a typical lamina within this laminate, the 
relationship between stresses and strains referred to x, s coordinate system may be 
expressed in the following matrix form 



Where e and e are the normal strains in x and s coordinate directions, y is 

X S 

the engineering shear strain along the contour of the cross section. Q„ are the 
stiffnesses which depend upon the fiber orientation and elastic moduli of the 
lamina. Equation (3) is derived in Ref. (12). 

Now, the layers are assumed to deform together so that the skin can be 
represented as a homogeneous anisotropic material represented by averaged 
material properties. Let N , N x and N x$ be the stress resultants per unit length in 
s-coordinate direction. Then, 
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Where 


X' 


— • 
<jT 

1 

A 12 

i 

N , 

> = 

A 12 

A 22 

A 26 

N « 


- A 16 

A 26 

A 66- 



(4) 



i,j=l,2, and 6 


(5) 


and N denotes the number of layers in the skin. N is negligible in unpressurized 
slender structures, so it is set to zero. Equation (4) is contracted to the following 
results 



A J1 " A 12 A 12 / A 22 
A 16 ' A 12 A 26 /A 22 


A 16 “ A 12 A 26 /A 22 

M 

A 66" A 26 A 26 /A 22. 

W 


( 6 ) 


It is convenient to have Equation (6) expressed in the foliowing form: 


r N x] 


' H ll 

i — 

«« 

CM 

X 

- 

► - 



l Y xsJ 


, H 21 

H 22j 



(7) 


where 

H 22 = *^ A 66 " A 26 A 26^ A 22^ 

H 21 = ‘ tA I6 “ A I2 A 26^ A 22 )H 22 

Hjj = - A 12 A|2/A22 + (A J6 - A|2 A 26^ A 22^ H 21 

(8a) 

The resultant loads acting on a given cross section are defined by the 


following integrals; 
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(N, M y , M z ) = J> N x 0, 2 , y)ds (8b) 

Where, N Is the tensile load along x-axis, M y is the bending moment about y-axis 
and M z is the bending moment about z-axis. 

The location of the origin of the y,z axes and their orientation are at 
present arbitrary. The ensuing analysis is greatly simplified if they are selected so 


that 



(y, z, yz)Hjjds 


( 0 , 0 , 0 ) 


(9) 


The resulting locus of origins of each section, then, defines the tension axis of the 
blade and the new y, z axes are the principal axes of bending. Now the stress 
resultants are evaluated with the aid of equations (1) and (7). 



( 10 ) 

(ID 

( 12 ) 


Since the cross sections undergo only a rigid body rotation and also because 

the axial displacement, u, is continuous along the contour of the cross section, the 

(13) 

rate of twist is evaluated from the following formula' ' 




»x 


and Equations (1) and (7) to obtain 



(13) 


Equations (2), 


= <u ,xj 

fH 21 ds - W ,xxi 

% 

H 2 jZ ds 


" V ,xxj 

h “* 

^H 22 ds)/2S 

(14) 


(10), (11), (12) and (14) are now recast ir. the following form to obtain 
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the overall engineering elastic law for the structure 


' u ,*' 


" s .i 

0 

0 

S l<f 


r N " 


- w ,« 

< 

►= 

0 

S 22 

0 

S 24 

i 

M y 

► 03) 

"\xx 


0 

0 

S 32 

S 30 


M z 


♦ f K 

V ' J 


> 

S 2k 

S 34 



M x 

V. > 



S is a symmetric flexibility matrix whose nonzero elements are given by 

s n = 

*1* * V (25! / H 2l ds) 

S 22 = 1 /^H|jZ 2 ds 

S 24 = S 22^ 2 ^/ H 21 z ds ^ 

s 33 “ 1//H n y 2 ds 
S 34 = S 33 /(2iT^H 22 y ds) 

= ^2^ H 22 ds + Jjj- S I^H 2 l ds 

+ S 2^/ H 21 zd5+ ^ S 2^f/ H 21 yds (l6) 

The development described so far is reported in Ref. (1). The constitutive 
equations (15) are to be supplemented with the following equilibrium equations 


N ,x + \ 


= 0 


M 


y»* 



(17) 


( 18 ) 
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z ,x 1 Q y 

(19) 

1 + q = 0 

z,x 

(20) 

V,x + <4 " 0 

(21) 

'x,x* m x * 0 

(22) 


to form the engineering theory for singie cell laminated beams under combined 

loading. Here, n represents the intensity of axial distributed loading, Q and Q 
x t y 

represent shear stress resultants, q and q represent the distributed loads and m 

£ y " 

represents the distributed torque. 

In order to proceed with the development of the new theory, it is necessary 
to have statically equilibrating stress resultants. The axial stress resultant can 
be represented in terms of overall stress resultants with the aid of liquations 05), 
(1) and (7) 

N x = fjN + f 2 M y + f 3 M 2 + f^M x (23) 

where 

fj = HjjSj, 

f 2 = z H 1 |S 22 

f 3 = y HjjS^ 

*4 = H 11 S 14 + z H 11 S 24 + y H 11 S 34 " H 21^ 2 ^ 

(24) 

The nonuniform shear flow due to bending and stretching can be derived from the 
following equilibrium equation 


N 


x,x 


+ N 


xs,s 


0 


(25) 


li 


and an additional requirement that the shear flow be equivalent to the applied 
torque. This result can be given as 


= A,x ds+ ^f ds K« d! < 26) 

With r denoting the perpendicular from the origin to the tangent at s. The axial 
stress term can be eliminated from above with the aid of Equation (25). 


where 


M 

N xs ■ - h N ,x - f 6 M y,x - '7 M x,x ‘ f S M x,x 


( i * /'«-M ds ' "g/ r ds /‘(i-l 0 ds ! i=5 ' 6 ’ 7 ’ 8 


(27) 

(28) 


Equations (23) and (27) provide stress resultants which are in static equilibrium. A 
comparison of Equations (2) and (27) reveals that the two expressions for the shear 
flow are different. The new shear flow is nonuniform along the contour of the 
cross section and this is due to bending and stretching loads. The nonuniform shear 
flow in turn affects the twist; its contribution can be evaluated from Equations (13), 
(7), (15) and (28). The result can be expressed as 

*,x = V +S 2A + S 34 M z + Vx 


f 9 N ,x + f 10 M y,x + f ll M z,x + £ 12 M x,x 


where 


(29) 


f l = - 7^ H 21 W ' U9 ' I0 ' H ’ 12 <30) 

fjQ and f|j will be zero if the shear resultants pass through the shear center. This 
does not, however, indicate that twist will not develop under bending and stretching 
related loads. The presence of coupling includes this possibility. 
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DEVELOPMENT OF THE NEW THEORY 
The theoretical developmental process is described in the form of a flow 
chart in Figure 4. The fundamental assumption that permits this development is 
that the transverse shear strain can be estimated from the statically equivalent 
stresses given by Equations (23) and (27). It is also assumed that the cross sections 
undergo a rigid body rotation. If n is the tangential displacement in the s- 
direction, then 

n = r 4> (31) 


The linear strain-displacement relation for the transverse shear strain is given by 


Y xs = u ,s + V 


(32) 


The following result emerges upon integration of u s 

u = U~zW -yV -2D<J) 

7* 7 X jX 


with the aid of Equation (31) 
+ / Y xs d3 U3) 


where 



(34) 


The first three terms in Equation (33) satisfy the engineering hypothesis that plane 
cross sections remain plane during deformation. The fourth term denotes the 
effect of restrained warping and the fifth term contains the effect of transverse 
shear deformation. With the aid of Equations (7), (10, (15), (27) and (29) the non- 
classical part of the displacement is expressed in terms of force resultants to 
arrive at 


u = U - zW >x - yV x + f, 3 N - f w M y ♦ fj 5 M z + * 16 M X 


+ f 17 N ,x + f 18 M y,x + f 19 M z,x + f 20 M x,x 


(35) 


where 
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fj3 = -2a,S 14 + / H 21 S ll ds 
f 14 = - 2w S 24 + y* H 21 S 22 2 ^ 

f l5 s + / H 2i S 22 yds 

fj£ = - 2u> + > y'^21^14 + ^21 Z ^24 + ^31 ^34 + ^22^ 2 ^ ds 

f 17 = - 2ti)f 9 -/H 22 f 5 ds 

f 18 = " 2ft) f 10 "/ H 22 f 6 ds 
f 19 = " 2u> f ll _ / H 22 f 7 ds 
f 20 = ~ 2a) f J2 " t /' H 22 f 8 ds 

The new displacement is now used in conjunction with Equations (7) and (26) 
to compute the axial stress resultant 


N. 


= H ll (U ,x 


zW ,xx - yV ,xx 


♦ J| 3 N x * f,«M yjX ♦ f,jM 


z,x 


+ f 16 M x,x 


+ f 17 N ,xx + f 18 M y,xx + f i9 M z,xx + f 20 M x,xx > 

M 

- H 21 ( * f 5 N ,x - f 6 M y,x " f 7 M z,x “ f 8 M x,x> (37) 


The engineering kinematic variables U, W, V must be related to the generalized 
force variables. This is accomplished by enforcing the definitions in Equation (8) 


<N,M f M z ) = U x /h u (I, z, y)ds 

■ W ,xx/ H n <1 ’ z ’ y) zds 

- V ,xx/ H ll U - z ’ y)yds 
+ N J [H n f i3 + H 21 f 5 )(l, z, y)ds 


14 


* M y, J 

^ H ll f l4 + H 21 f 6^ 1 ' 

z, y)ds 


^ H ll f 15 + H 21 f 7^» z > y * ds 

* M x,yJ 

r (Hji f i 6 + H 21 f 8 )(l, 

z, y)ds 

♦%/ 

Huf i ?(!» y, z)ds 


* M y,xx. 

/H,if 18 (l’ y, z)ds 


+ M Z,XX 

|'H 11 f jg (l, y, z)ds 


+ M x,xx, 

f ^u f 20 (1, y ' z)ds 


+ M x /(-H 2 ,/2 flXl, y,z)ds 



These equations can be cast in a more familiar form by solving them for the 
engineering kinematic variables. 
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It 


N \ 
»xy 1 

5 H 

5 12 

S 13 

S 14 


It 

M 

It 

II 


^y,xx( 

S 12 

S 22 

S 23 

S 24 

* 

ft 

ft 

tl 

ft 


M z,xxl 

5 13 

S 23 

S 33 

s 34 


ft 

It 

ft 

II 


( M x,xx/ 

S 14 

— 

S 24 

s 34 

S 44 



(39) 


Equations (39) describe the stretching and bending of the beam. They replace the 
classical Equations (15). The additional flexibility matrices Cs!)and Cs"j depend only 
on the elastic modulii of the material and dimensions of the cross section. The 
expressions for these parameters can easily be deduced. For example 

s'n ■ <4 H jjdsV^Hjjfjj + H 21 f 5 )ds, etc. (40) 

The axial stress resultant may be expressed in terms of generalized force 
variables with the aid of Equations (39) and (37) 

N x r fjN + f 2 M y + f 3 M z + f^M x 

+ f 21 N ,x + *22 M y,x + f 23 M z,x + f 24 M x,x 
+ f 25 N ,xx + f 26 M y,xx + f 27 M z,xx + f 28 M x,xx 

(41) 


f 21 = H 11 (S 11 + zS 12 + yS 13 + f 13 } + H 21 f 5 
f 22 ~ H 11^ S 12 + zS 22 + yS 23 + + H 2i f 6 

f 23 = H 11^ S 13 + zS 23 + yS 33 + f 15^ + H 21 f 7 

f 24 = H 11 (S 14 + zS 24 + yS 34 + W + H 21 f 8 

f 25 = H ll (S "ll + zS "l2 + yS 'l3 + f 17> 


where 


16 


26 

■ h ii (s 

12 

+ zS 

22 

♦ yS 

23 

+ f 18 ) 


1 

1 

If 

II 

+ f i9 ) 

27 

= H|,<S 

13 

+ zS 

23 

+ yS 

33 


♦ 1 

II 

ft 

+ f 20 ) 

28 

= H n (S 

14 

+ zS 

24 

+ yS 

34 


(42) 


The corrections to the axial stress resultant include the effects of bending 
related warping and torsion-related restrained warping. The new shear flow can 
now be estimated substituting this stress in Equation (26). 


N 


xs 


M. 


2ft 


" f 5 N ,x " f 6 M y,x " *7 M z,x ' f 8 M x,x 


' f 29 N ,xx ’ f 30 M y,xx ’ f 31 M z,xx " f 32^x,xx 


' f 33 N ,xxx ' f 34 M y,xxx ' f 35 M z,xxx ’ f 36 M x,xxx 

(43) 

where 

f i * /*(i-8) ds - ds ‘= 29 ' 30 - ' ' ■ 36 

The rate of twist is re-evaluated with the aid of Equations 03), (7), (35) and (41) 


♦,x = S 14 N +S 24 M y + S 34 M z + 5 44 M x 

+ *9 N ,x + *lO M y,x + f li M z,x + f 12 M x,x 

+ f 37 N ,xx + f 38 M y,xx + f 39 M z,xx + f 40 M x,xx 
+ f 41 N ,xxx + f 42 M y,xxx + f 43 M z,xxx + f 44 M x,xxx 


where 
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*37 2 jj ^ ds|n 21 (S u + zS 12 + yS 13 + i [? ) - H 21 f 29 j 

*38 = / ds | H 21 (S 12 + zS 22 + y S 23 + f l8 ) “ H 21 f 3o| 

*39 = 7S f ds | H 21 (S 13 + zS 23 + yS 33 + f 19 J " H 21 f 31 1 

f ^° = $ ds | H 2i^ S 14 + zS 24 + yS 34 + W " H 2l f 32j 

f i = j> ds |- H 2i f (i-8)| * 1 a W » k2 ’ 43 > 44 <46) 


The governing equations for the response of single cell laminated rotor 

blades can now be summarized. Overall equations of equilibrium remain same as 

those in (17) - (22), Similarly, the generalized kinematic variables and force 

variables also remain the same. The kinematic variables are U x , -W ^ -V ^ and 

$ . The force variables are N, M , M and M . Equations that relate these two 
,x y z x 

sets of variables are given by (39) and (45). These overall constitutive equations 
contain the effects due to transverse shear, torsional shear and restrained warping. 
New distributions for the axial displacement, axial stress and shear flow are 
contained in Equations (35), (41) and (44). 

APPLICATION 

The following section provides an example which illustrates the relative 
significance of nonclassicai influences on overall response predictions and stress 
estimates. The potential for aeroelastic tailoring is also investigated to a limited 
extent. For this purpose, the rotor blade of Boeing’s CH-47 helicopter is used as a 


18 


guide in sizing. The blade is modeled as a single cell beam with a uniformly 
distributed load, q z , with q applied at the quarter-chord point. The end where the 
blade is joined to the hub is treated as a clamped end, with the other end being 
free. The dimensions of the box are shown in Figure 5. 

In order to study the influence of material properties, two different 
material systems, graphite/epoxy and glass/epoxy, are selected. The following 
material properties are assumed for each ply: 

Graphite/Epoxy; 


Ej = 25 x 10 6 psi, 
G 23 = 0.2 x 10 6 psi, 


\> 


12 



0.25, 


E 2 = E 3 = 1 x J0 6 psi 
G 12 = G 13 = 0.5 x 10 6 psi 
v 23 = 0.25 


Glass/Epoxy : 

E, = 7.3 x I0 6 psi, E 2 = E 3 = 2.1 x 10 6 psi 

G |2 = Gj 3 = 0.88 x 10 6 psi 
v 12 =v 13 = 0.275, v 23 = 0.25 


The E's are Young's modulii, the G's are shear modulii and v's are Poisson ratios. 
The subscript "1" denotes a direction along the fibers, "2" and "3" denote directions 
normal to the fiber. All plies are identical with a thickness of 0.0085 inches. 

Two different layups are considered for the box skin. These are (0j 2 , 45 j 2 ) 
and (0^, 45g). The first one is more effective in torsion and the second one in 
bending. Furthermore, these plies are arranged with respect to z = 0 plane in such 
a way that there exists three distinctly cifferent couplings between elastic modes of 
deformation. Because of the nature of the applied loading and skin layup, there is 
no coupling between stretching and twisting or between bending about z-axis and 


19 


twisting. The only coupling is between bending about the y-axis and torsion (flap- 
wise bending). This coupling is controlled by S 2(f . As shown in Figure 6, three 
different layups are possible depending upon whether 45° plies are placed 
symmetrically or antisymmetrically with respect to z = 0 plane. 

The first configuration is a balanced design. Since the 45° plies are placed 
symmetrically, is zero, which means that there is no coupling between bending 
about y-axis and torsion. The other two designs are unbalanced. Depending on the 
angle of orientation with respect to the load, these designs produce either enhanced 
(positive) coupling or reduced (negative) coupling. As shown in the Tables 1-3, there 
are now twelve different designs depending upon the material system, relative ratio 
of 0° and 45° plies and relative orientation of the 45° plies in upper and lower parts 
of the box and placement of the load. 

Under the applied uniformly distributed line load (q = q), all these 
structures undergo twisting as well as bending about y-axis. Since this is a 
statically determinate problem, the flapwise bending moment and torque are 
estimated as 

M y = -§(x-L) 2 

M x = -q e(x - L) (47) 

e is eccentricity of applied loading which in this case is 5 inches, The applicable 
boundary conditions are given by the following: 

At x = 0: W = 0, W v = 0, <|> = 0 (48) 

Under these conditions, equilibrium equations (18), (20), (22), the overall 
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constitutive equations for -W from (39) and the Equation (45) for twist are 
solved. 

The results for the transverse displacement and the twist at the free end 

are presented in Tables 1 and 2, respectively. These parameters are given in terms 

of unit applied load, q. Similarly, the maximum values for the stress resultants 

and N at the root are also presented per unit applied load in Figure 3. The 
xy 

maximum value of the axial stress occurs at the center of the web-like vertical 
portion of the structure. 

The present theory results for W and 4» are subdivided into three groups 
as shown in Tables 1 and 2 to illustrate the improvements due to nonclassical 
effects and coupling. In the simplest engineering approach, the coupling effects 
are usually neglected. The material is assumed to be orthotropic, effectively 
represented by equivalent bending and torsional rigidities. These results are 
denoted "classical" in Tables 1 and 2. It can be concluded from this study that the 
nonclassical effects influence deflection more than the coupling parameter and the 
improvements from both coupling and nonclassical effects are relatively small 
compared to the classical result. On the other hand, the coupling accounts for a 
major portion of the twist. If it is not incorporated, major discrepancies arise in 
the twist predictions. This is because the classical approach cannot delineate 
between the balanced and unbalanced designs. The nonclassical Influences on twist 
predictions are negligible. With respect to the transverse displacement prediction, 
it may be seen that the nonclassical influences are more pronounced for 
graphite/epoxy material system and (0 ^ 45g) layup. 

Table 3 shows predictions for the maximum values of the stress resultants. 

Present theory results are evaluated from Equations (41) and (43). For the sake of 

2 

comparison, stresses from Mansfield's theory are also presented. These are not 
available in Ref. (2), but are evaluated using Equations (1), (7) and (15). Unlike the 
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"classical" results for twist and deflection, Mansfield stresses include the bending- 
torsion coupling effect. Therefore, the apparent differences between the present 
theory predictions and those from the Mansfield theory are entirely due to 
nonclasslcal influences. Failure to account for these effects leads to an under 
estimation of the stresses. This could be as much as 4% for tensile stresses and 9% 
for the compressive stresses. 

The structure considered here is very slender. The theory is also developed 
under the premises of a thin walled cylinder. This explains why the nonclasslcal 
influences are so small in this example. Experience^® shows, for beams and 
plates made of composite materials, the nondassical influences become important 
enough to be included as the structure becomes "thicker" on a relative basis. Since 
the kinematic procedure used here is the same as the one used in Ref. (7)-0O), the 
present theory can be expected to yield valid results even for very thick-walled 
rotor blades. 


SUMMARY AND CONCLUSIONS 

A theory has been developed for composite rotor blades that can be 
modeled as dosed cell beams. The theory accounts for nonclasslcal influences 
related to restrained warping, transverse shear stain and bending and stretching- 
related warping. A primary influence that is included is elastic coupling that 
results from arbitrary composite ply layup in the walls of the cell. To illustrate the 
role of elastic coupling and nondassical influences, a rotor blade-type model has 
been created. This model is analyzed as a single cell cantilever beam under 
uniformly distributed load. The following conclusions are drawn from this example: 
1. Nondassical influences and elastic coupling have a negligible effect 
on transverse deflection. 

Twist is controlled by the elastic coupling. Nondassical influences 
have a negligible impact on twist. 


2. 
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3. The axial stress predictions are improved as a result of accounting for 
the nondassical Influences. The shear flow is affected only slightly 
by these nondassical influences. 

4. The potential for aeroelastic tailoring has been illustrated by 
choosing three designs with different layups. The first one is a 
balanced design with zero coupling. The second and third are 
unbalanced designs; one of these two produces enhanced coupling and 
the other reduced coupling. For all of these three designs, the weight 
is the same and transverse deflection and maximum stresses 
generated within are also approximately the same, yet the induced 
twists are different. 

5. The structure considered in this example is very slender. This is the 
reason why the nondassical influences appear so small in this 
example. By the experience garnered by the authors with laminated 
thick plates and beams, it is believed that these influences become 
more pronounced for thick walled rotor blades and the present theory 
is capable of predicting them. 

RECOMMENDATIONS 

It appears that the nondassical influences are of little importance in the 

model. As the model approximates practical dimensions of main rotor blades, it is 

likely that the slenderness, together with dosed cell construction, are responsible 

2 

for this. As a practical matter, the dassical theory of Mansfield and Sobey should 
be adequate in most cases and should be used. 

If the blade mder consideration is less "slender" than the example or if 
moisture-temperature effects significantly alter the stiffness properties of the 
composite material, then the new theory developed herein can be used. Another 
situation, the determination of higher than first few vibration modes, also requires 
the use of the present, more complicated theory. 
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FIGURE 6. SKIN LAYUPS 
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Table 1. Transverse Displacement at the Free End 





























































Table 2. Twist at the Free End 
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